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Abstract 

We discuss the behavior of charmonia and bottomonia correlators and spectral func- 
tions above the deconfinement temperature and determine melting temperatures for 
different mesonic states. 
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1 Introduction 

The study of heavy quarkonia at finite temperature is an important way to 
probe hot and dense strongly interacting matter. General considerations sug- 
gest that quarkonia will melt at temperatures somewhat higher than the de- 
confinement temperature as a result of modification of interquark potential 
due to color screening. Matsui and Satz [1] proposed that color screening at 
high temperature will lead to quarkonium suppression, which can be used as 
a signal of QGP formation. Studies of modification of quarkonium potential 
in terms of spectral functions have started to appear only recently [2,3,4,5]. 
However, the systematic errors were not well controlled and it is not clear 
which features of the charmonium spectral functions are physical and which 
are lattice artifacts of the MEM analysis and/or finite lattice spacing. 

This contribution reports updated results of our study of charmonium [6,7] 
and bottomonium [8] spectral functions and correlators at finite temperature. 
The calculations are performed in quenched QCD using anisotropic lattices, 
the Wilson gauge action and the Fermilab fermion action [11]. Further details 
of the lattice action and the parameters of the simulations can be found in 
[6]. Detailed results of the study and description of the Maximum Entropy 
Method (MEM) will be presented in our upcoming paper [9]. 
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2 Meson Correlators and Spectral Functions 



We look at correlators of point meson operator JuitjX) = q(t,x)THq(t,x), 
where = 1, 75, 7^, 757/^, 7^7i/ fixes the mesonic channel to (pseudo) scalar, 
(pseudo)vector and tensor channels. On the lattice one computes the Euclidean 
propagator Gh(t,'p) = J (Pxe^^'^ {T^-Jh^t, x)Jh{0,0)) , which is related to the 
spectral function through the integral representation 



G{t,p) = J diuaK{uj,T), K{uj,t) 



cosh(cj(r - 1/2T)) 
sinh(cj/2T) ' 



Spectral functions can be reconstructed using Bayesian techniques, such as 
the MEM [10]. Instead of using Bryan's algorithm we use a new algorithm, 
which will be described in details in [9]. 

First we check the MEM ability to reconstruct a model spectral function of 
finite width. For a model spectral function we choose the Breit-Wigner model 
spectral function, which (in physical units) is given by 

, , 1 7 3 
a{uj) = ri— — — J—*—- 2 



The covariance matrix from an actual lattice simulation of charmonium at 
zero temperature (/5 = 6.5, ^ = A, Nt = 160) is used to introduce the noise. 
In Fig. 1 we show the model and reconstructed (with mioj) = O.OSSiu;^ and 
m{u) = 1 default models) spectral functions for the pseudoscalar channel. We 
see that the MEM can successfully reconstruct the model spectral function. 
However, it is gradually loosing its ability as we increase the width and at 
7 = 2 GeV is incapable of spectral function reconstruction. 
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Fig. 1. Reconstruction of the 7 = 0.1 and 0.4 GeV Breit-Wigner model spectral 
functions. 
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3 Spectral functions and correlators at finite temperature 



With increasing temperature it becomes more and more difficult to recon- 
struct the spectral functions, therefore it is useful to study the temperature 
dependence of quarkonia correlators ffist. To separate out the trivial tem- 
perature dependence due to the integration kernel, following Ref. [5] at each 
temperature we calculate the so-called reconstructed correlator 

OO 

Greconir, ^) = / duja{uj, T = 0)K{t, lu, T) (3) 





The ratio of the original and the reconstructed correlator should be close to 
one, G{t,T) /Grecon{'r,T) ~ 1, if there is no temperature dependence in the 
spectral function. In Fig. 2 we show the ratio of the correlators for charmo- 
nium in pseudoscalar and scalar channels. It is evident that there is little 
modification of the correlator in the pseodoscalar channel for temperatures 
up to 1.5Tc. For the scalar channel, however, the change in the correlator is 
noticeable at temperatures as low as 1.16Tc. The results for (3 = 6.1 and 6.5 
are very similar, thus there is little dependence on the lattice spacing. In Fig. 
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Fig. 2. The ratio G{t,T) / Grecon{T-,T) of charmonium for pseudoscalar (left) and 
scalar (right) channels at aj"^ = 8.18 and 14.11GeV at different temperatures. 

3 we present charmonium and bottomonium spectral functions for the scalar 
channel at different temperatures. From these figures the dissolution of the 
ground states at temperatures above deconfinement is clear. For charmonium 
we tried two different default models m{u) = 0.038ci;^ and m(a;) = 1 and the 
results are qualitatively the same, which is another check on systematics. 



4 Summary 

From our studies of correlators and spectral functions we can conclude that, 
the IS {rjc, J/'ip) charmonium states exist as a resonance in the deconfined 
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Fig. 3. Charmonia (left) and bottomonia (right) spectral functions in scalar channel 
for different temperatures. 

phase at T ~ 1.5Tc. On the other hand, the IP {Xco, Xa) charmonium states 
dissolve at l.lTr. Bottomonium states show similar behavior. 
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